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ABSTRACT 



We present a kinematic analysis of the globular cluster (GC) system in the 
U ■ giant elliptical galaxy (gE) NGC 4636 in the Virgo cluster. Using the photometric 

and spectroscopic database of 238 GCs (108 blue GCs and 130 red GCs) at the 
galactocentric radius 0'.39 < R < 15'. 43, we have investigated the kinematics of 
the GC system. The NGC 4636 GC system shows weak overall rotation, which is 
dominated by the red GCs. However, both the blue GCs and red GCs show some 
rotation in the inner region at R < 4'. 3 (= 2.9-R e fj = 18.5 kpc). The velocity 
dispersion for all the GCs is derived to be a p = 225+g 2 km s _1 . The velocity 
dispersion for the blue GCs (cr p = 251^ km s _1 ) is slightly larger than that for 
CN . the red GCs (a p = 205^ km s _1 ). The velocity dispersions for the blue GCs 

about the mean velocity and about the best fit rotation curve have a significant 
. variation depending on the galactocentric radius. Comparison of observed stellar 

and GC velocity dispersion profiles with the velocity dispersion profiles calculated 
from the stellar mass profile shows that the mass-to-light ratio should increase 
as the galactocentric distance increases, indicating the existence of an extended 
dark matter halo. From the comparison of the observed GC velocity dispersion 
profiles and the velocity dispersion profiles calculated for the X-ray mass profiles 
in the literature, we find that the orbit of the GC system is tangential, and that 
the orbit of the red GCs is slightly more tangential than that of the blue GCs. 
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We compare the GC kinematics of NGC 4636 with those of other six gEs, finding 
that the kinematic properties of the GCs are diverse among gEs. We find several 
correlations between the kinematics of the GCs and the global parameters of 
their host galaxies. We discuss the implication of the results for the formation 
models of the GC system in gEs, and suggest a mixture scenario for the origin 
of the GCs in gEs. 

Subject headings: galaxies: clusters: general — galaxies: individual (NGC 4636) 
- galaxies: kinematics and dynamics — galaxies: star clusters 



Introduction 



Globular clusters (GCs) are an excellent tracer to probe the gravitational potential 
of their host galaxies. Kinematics of the GC system in a galaxy is determined by the 
gravitational potential of their host galaxies and the GCs contain a fossil record of the 
dynamical evolution since the formation of their host galaxy. From the kinematic study of 
the GC system we can estimate the global mass distribution of their host galaxy including 
the dark matter, or can derive information on the orbital properties of the GCs if a mass 
distribution of their host galaxy is known in prior (e.g., X-ray observation of the hot gas). 
GCs are useful especially for the study of the outer region of a galaxy where the gravitational 
potential is mainly dominated by the dark matter. They are particularly efficient for the 
study of nearby giant elliptical galaxies (gEs) where thousands of GCs are often found. 



Recently iHwang et al.l (120081 ) presented a good summary of the results on the kinematics 
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NGC 1407, in a nearby galaxy group and provided updated kinematic analysis of the GC 
systems in gEs using their consistent method. One surprising finding from these studies is 
that the kinematics of the GC systems in these gEs is diverse, showing a large difference 
in velocity dispersion, rotation, and their radial variation among the gEs. This is in stark 
contrast to the photometric prope rty of the GCs that is more or less similar among the gEs 



(1Ledl2003l : iBrodie fc Strader! 120061 ). 
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We have been carrying a project to investigate the spectroscopic properties of the GCs 
in nearby galaxie s. Our study on t he ki nematics of the GC system of M60, a gE in Virgo, 
was presented in ILee et al.l (I2008af) an d iHwang et al.l (120081 ). and another on the M31 GC 
system was given in ILee et al.l (l2008d ). Recently we presente d the measurement of radial 
velocities for the GCs in NGC 4636 in the companion paper (jPark et al.l l2009al . hereafter 
Paper I), and we present a detailed kinematic analysis of these data in this paper. 

NGC 4636 is an E/SO galaxy in Virgo. It is located 10°.8 (2.8 Mpc at the distance of 
NGC 4636) south east from the Virgo cent er, and is conside r ed to be a major mem ber of 
a small group falling into the Virgo center (jNolthenius Ill993l ; iTrinchieri et al.lll994l ). NGC 
4636 is relatively less luminous [My = —21.7 mag) among the gEs in Virgo. It shows several 
interesting structures in the central region: jets in the radio images, arm- like struc tures, 



bubbles, and a cavity in the X-ray images (jO'Su 



livan et al 



Temi et al. 



20051: iBaldi et al.l 120091 ). and 



2003h . It is notable that the 



dust emission in the far- infrared (100/xm) image 
faint isophotes in the outer region of NGC 4636 are much flatter (E4) th an those in th e 
inner region (E0), indicating that this galaxy may be in the transition to SO (ISandagelll96ll ). 
This fact may imply the presence of large-scale angular momentum associated with recent 
m ergers. We adopted a distance to NGC 4636, 14.7 Mpc [(m — M) = 30.83±0.13], as given 
by lTonry et al.l (120011 ) based on the surface brightness fluctuation method. One arcmin in the 
sky at this distance corresponds to 4.26 kpc. Effective radius, ellipticity, and position angle 

kpc, e e fr = 0.256, and PA c ff = 148 deg, respectively 
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This paper is composed as follows. Section [2] gives a brief description of the data used 
in this analysis, and the kinematic properties of the GC system in NGC 4636 are derived in 
§|3j In §H we compare kinematic properties of the NGC 4636 GC system to those of other 
gEs, and discuss the implication of the results regarding the GC formation models. Primary 
results are summarized in the final section. 



2. Data 

We used the spectroscopic data of the GCs given in Paper I, which describes the details 
of the spectroscopic observation, data reduction, and the data set. Here we only give a 
brief summary of the data set of NGC 4636 GCs. We selected GC candidates in deep, 
wide-field W ashing ton C and T% images (15.8' x 15.8') obtained at KPNO 4m telescope 



( IPark et al.ll2009bl ). Spectroscopic observations of these targets were made using the Multi- 
Object Spectroscopy (MOS) mode of Faint Object Camera and Spectrograph (FOCAS) on 
the SUBARU 8.2m Telescope. 
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We determined the radial velocities for 122 objects including 105 GCs, 11 foreground 
stars, 2 background galaxies, 3 probable intracluster GCs in Virgo, and the nucleus of NGC 
4636, by cross-correlating the target spectra with five Galactic GCs as templates. Paper I 
presented a master catalog of the radial ve locities for the GCs in N GC 4636 combining their 
data w ith the data for 1 74 G Cs given by ISchuberth et al.l (120061 ). The radial velocities of 
GCs in lSchuberth et al.l (120061 ) were transformed into our velocity system using equation (1) 
in Paper I, and the transformed velocities were used for further analysis. 

Of the entire spectroscopic sample of GC candidates we selected 238 member GCs of 
NGC 4636 in Paper I using radial velocities (300 < v p < 1600 km s _1 ) and (C — Ti) colors 
(0.9 < C - Ti < 2.1). There are 108 blue GCs (0.9 < (C - Ti) < 1.55) and 130 red GCs 
(1.55 < (C — Ti) < 2.1) in the total sample. In Figure [fl we show the spatial distribution of 
the objects in NGC 4636 with measured velocities. It shows that the spatial segregation of 
high-velocity GCs (open symbols) and low-velocity GCs (filled symbols) is not clearly seen, 
which indicates that the rotation of the GC system is, if any, weak (to be discussed in §3. 1[) . 



Results 



We have investigated the kinematic properties of the GC system using the master catalog 
of 238 GCs in NGC 4636 (Paper I): the rotation amplitude, the position angle of the rotation 
axis, the mean line-of-s ight velocity, the proje cted velocity dispersion, and the velocity ellip- 
soid. ICote et al.l (120011 ) and ICote et al.l (120031 ) presented detailed analysis of the kinematics 
of the GC systems in M87 and M49, respectively. We adopted the anal ysis method used b y 
these studies for the following analysis, as done for the M60 GC system ( iHwang et al.ll2008l ). 



3.1. Rotation of the GC System 

First we derive the rotational property of the GC system in NGC 4636. We made 
assumptions to derive the intrinsic rotational velocity field from the radial velocity data: 
(a) the GC system is spherically symmetric with an intrinsic angular velocity field stratified 
on spheres, and (b) that the GC rotation axis lies in the plane of the sky. With these 
assumptions, it is expected that radial velocities of GCs depend, if they follow any overall 
rotation, sinusoidally on the azimuthal angles. For the NGC 4636 GC system, the assumption 
of spherical symmetry of the GC system is reasonable due to the modest projected ellipticity 
(effective ellipticity, e cfT = 0.26). 

Thus we can determine the amplitude and axis of the rotation, fitting the observed 
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line-of-sight velocities (v p ) of the GCs with the function, 



v p (Q) 



J sys 



+ (SIR) sin(6 - e ) , 



where G is the projected position angle of a GC relative to the galaxy center measured from 
north to east, Go is the projected position angle of the rotation axis of the GC system, 
R is the projected galactocentric distance, (QR) is the rotation amplitude, and t> sys is the 
systemic velocity of the GC system. 

Figure [2] displays the radial velocities of GCs with measured uncertainties as a function 
of position angle for all 238 GCs (top), 108 blue GCs (middle), and 130 red GCs (bottom). 
We overlaid the best fit rotation curve of equation ([I]) for each sample. For fitting the data 
we used an error-weighted, nonlinear fit of equation (TjQ) with v sys as a fixed value of the 
velocity of the NGC 4636 nucleus (v ga \ = 928 ± 45 km s _1 ) rather t han as a fre e param eter 
for a better fitting. We estimated, using the biweight location of iBeers et al.l (jl990t ). the 
systemic velocity of the GC system to be t> sys = 949t}g km s" 1 (see §3.21) . which agrees with 
the velocity of the NGC 4636 nucleus. 



The values of the rotation amplitudes (QR) we derived are 37^ km s 1 for all the GCs, 
27tf 4 km s" 1 for the blue GCs, and 68^ km s" 1 for the red GCs. Thus, all the GCs and 
the blue GCs show little rotation, while the red GCs show a marginal hint of rotation. Our 



result s are consistent with those based on a smaller sample of 174 GCs in ISchuberth et al. 



(|2006j) who derived 28 ± 18 km s" 1 for all the GCs, 11 ± 27 km s^ 1 for the blue GCs, and 
87 ± 18 km s" 1 for the red GCs. 

The orientation of the rotation axis (Go) is estimated to be 174°ljg for all the GCs, 
0°li44 for the blue GCs, and 178°^34 for the red GCs. The orientations of rotation axes 
for all subsamples appear to be similar, and they are closer to the photometric major axis, 
showing that the GCs ro tate around the major a xis (Q p hot = 148°). This result is in contrast 
with the result given by ISchuberth et al.l (120061 ) who derived 63° ± 43 for all the GCs, and 
60° ± 22 for the red GC s, showing that t he GC s rotate around the minor axis. If we use 
only the data g i ven b y ISchuberth et al.l (120061 ). we get similar values to those given by 
Schuberth et al.l (120061 ). Therefore, the difference in the orientation between the two studies 
is due to the difference in the sample. 

In Figure [3], we present the rotation of the GC system for the samples of different radial 
bins in order to investigate any radial variation of rotational properties. We divided the GCs 
into two groups: the GCs in the inner region at 23" < R < 260" and those in the outer region 
at 260" < R < 926". Each of the sub regions includes a similar number of GCs (120 and 
118, respectively). We applied the same fitting procedure as used for the results in Figure [2] 
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for all the GCs, the blue GCs, and the red GCs in each radial bin. 

Interestingly both the blue GCs and red GCs show some rotation around the major axis 
in the inner region (119lg2 km s_1 ' an d 77^ km s -1 , respectively), while both subsamples 
show little rotation in the outer region. The rotation of the red GCs seen for the entire 
radial range is found to be mainly due to the red GCs in the inner region. Little rotation is 
seen for all the GCs even in the inner region, although both the blue GCs and red GCs show 
measurable rotation in the inner region. This is because the blue GCs and red GCs rotate 
in the opposite direction to each other, canceling the rotational effect for all the GCs. 



3.2. Velocity Dispersion of the GC System 



We summarize the kinematics of the NGC 4636 GC system derived in this study in 
Table HJ Several kinematic parameters for all the GCs, the blue GCs, and the red GCs 
are presented for the entire region, the inner region, and the outer region. The column (1) 
defines the range of the projected radial distance from the center of NGC 4636 for each 
region in arcsec, and the column (2) gives the mean value of the radial distance in arcsec. 
The number of GCs in each region is shown in the column (3). T he columns ( 4) an d (5) 
represent the mean line-of-sight velocity (the biweight location of iBeers et al.l 119901) an d 
the velocity dispersion about this mean velocity (the biweight scale of IBeers et al.lll990l ). 
respectively. The position angle of the rotation axis and the rotation amplitude estimated 
using equation ([1]) in each region are given in the columns (6) and (7), respectively. The 
column (8) gives the velocity dispersion about the best fit rotation curve (cr p , r )- The column 
(9) gives the absolute value of the ratio of the rotation amplitude to the velocity dispersion 
about the best fit rotation curve. The uncertainties of these values represent 68% (la) 
confidence int ervals that are det ermined from the numerical bootstrap procedure following 
the method of ICote et all (boOlh . 



The velocity dispersion for all 238 GCs is estimated t o be 225lq 2 km s 1 . Th is result is 



consistent, within the uncertainty, with the value given by lSchuberth et al.l (120061 ) who used 
174 GCs, 203 ± 11 km s -1 . In addition, it is found that the velocity dispersion about the 
mean velocity of the GC system for the blue GCs (er p = 2511*2 km s _1 ) is 46 km s _1 larger 
than that for the red GCs (a p = 2051J 1 . km s" 1 ). We found a similar result for the velocity 



dispersion about the best fit rotation curve: a 



p.i 



252+^ km s" 



for the blue GCs, and 



a p>r = 203j^jg km s 1 for the red GCs. Our result for the red GCs is similar to the value 



given by lSchuberth et al.l (120061 ). a p = 199±17 km s 1 , within the uncertain ty. However, our 



value for the blue GC is approximately 50 km s 1 larger than that given by lSchuberth et al. 



(l2006h . o p = 202 ± 15 km s" 1 . 
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Figure H] displays the radial velocities of GCs with measured uncertainties against pro- 
jected galactocentric distances. The mean radial velocities in four radial bins are overlaid by 
squares with long horizontal error bar. The velocity dispersion about the mean velocity in 
each bin is also represented by a vertical error bar. The mean velocities of all samples agree 
well with the systemic velocity of NGC 4636. 

To investigate the radial variation of velocity dispersion in detail, we present a smoothed 
radial profile of velocity dispersion about the mean radial velocity (filled symbols) and about 
the best fit rotation curve (open symbols) in Figure [5J We calculate the velocity dispersion 
of the GCs lying within a radial bin with fixed width, AR = 120" ~ 8.52 kpc as increasing 
the bin center by a fixed step width, 5R = 10" ~ 0.71 kpc. We set the radial width and the 
step width so that the number of GCs per bin exceeds 10, and the calculation stops when 
the number of GCs in a bin is smaller than 10. 

The velocity dispersion for all the GCs varies more or less (about 50 km s" 1 ) around 
the mean value for the range of radius. However, the velocity dispersion for the blue GCs 
shows an abrupt decrease from « 300 km s~ x to ~ 150 km s _1 at R « 3', and then keeps 
increasing with increasing radius to ~ 300 km s _1 at R ~ T. On the other hand, the velocity 
dispersion for the red GCs stays almost constant at « 200 km s _1 from the center out to 
R ~ 6', shows a slight decrease of ~ 50 km s" 1 for 6' < R < 7', and increases slightly at 
R « 7'. The velocity dispersions about the best fit rotation curves of all three samples are 
not much different from those about the mean radial velocities (also see Tabled]). 



3.3. Velocity Anisotropy of the Globular Cluster System 

Assuming that the NGC 4636 GC system is spherically symmetric, we can apply the 
Jeans equation in the absence of rotation to the dynamical analysis of the GC system. The 
spherical Jeans equation is represented by 



d / \ 2 / \ i 2/? c i(r) 2 GM tot (r) 

— n cl (r)a r {r) H n d {r)a r {r) = -n d (r) 

dr r r z 



where r is a three dimensional radial distance from the galactic center, n c \(r) is a three dimen- 
sional density profile of the GC system, a r (r) is a radial component of velocity dispersion 
(5 c \{r) = 1 — V(j(r)/crf.(r) is a velocity anisotropy, G is the gravitational c onstant, and M tnt {r) 



is a t otal gravitating mass contained within a sphere of radius r (e.g., iBinney &: Tremaine 



19871 ). (Te{r) is a tangential component of velocity dispersion that is equal to an azimuthal 



component of the velocity dispersion, crAr), in the absence of rotation of the GC system. 
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There are two approaches to the dynamical analysis of the GC system: a) to derive 
the mass profiles of the host galaxy for given orbits of the GCs, or b) to determine the 
orbital properties of the GCs for a given mass profile of the host galaxy. Several studies 
on the dynamics of the GC system have focused on determining the gravitational mass, 
Mt t(r), using the Jeans equation by assuming a simple isot r opic orbit with c ](r) = (e.g., 
Cohen & Rvzhov!ll997l : iMinniti et aDll998l : IZepf et al.llioool : ISchuberth et~aDl2006h . 



However, with an aid of an inde pendent determination of the mass profile of an el 



liptical galaxy using X-ray data (e.g., iBrighenti fc Mathewd 119971 ; iMatsushita et al.l Il998 



Loewenstein fc Mushotzkv 



2003 



for NGC 463 6 ), the velocity an i sotropy itself can be inves - 



tigated (e.g., iRomanowsky fc Kochanekl l200ll ; ICote et al 



2001 



Cote et al 



op; 

gOpa iHwang et al.l l2008h . 



EST the M49 GC system 



Fo llowing; the analys is of the M87 GC system by 

by ICote et al.l (120031 ) . and the M60 GC system by IHwang et al.l (120081 ). we first derive the 
three dimensional density profile of the GC system, n c \{r) and the total mass profile, M tot (r). 
Comparing the velocity dispersion profile (VDP) calculated from the Jeans equation with 
the observed VDP, a p (R), we determine the velocity anisotropy of the NGC 4636 GC system. 



3.3.1. Density Profiles for the GC system 



Here we derive the three dimensional density pro files from the surfac e number density 
of NGC 4636 GCs for two models: the NFW profile (jNavarro et al.lll997l ) and the Dehnen 
profile (iDehnenl I1993T). W e used the surface number density profiles of NGC 4636 GCs 



Park et al.l J2009bh . They derived the surface density profile of NGC 4636 GCs by 



given m 

combining the HST/WFPC2 archive data for the inner region at R < 1.5', and the KPNO 
data for the outer region at R > 1.5'. They a dopted the background levels of the mean 
surface number density from lDirsch et al.l (120051 ) at R > 13.5': 2.33 ±0.07 per square arcmin 
for all the GCs, 1.88 ± 0.06 per square arcmin for the blue GCs, and 0.76 ± 0.04 per square 
arcmin for the red GCs. Then they subtracted these background values from the original 
number counts to produce the radial profiles of the net surface number density of GCs. 
Since they selected GCs that are brighter than T\ m 23.0 mag, it is needed to correct the 
surface number density profile in order to account for the uncounte d GCs due to the lim iting; 
magnitude. To calculate the correction factor, the equation (11) in iMcLaughlinl (119991 ) with 
Vum.x = oo was used on the assumption that the GC luminosity function of NG C 4636 has a 
Gaus sian shape with a peak at 7\ m 23.31 mag and a dispersion a = 1.03 mag dDirsch et al. 



20051). It is found that the surface number density of the bright NGC 4636 GCs in 
( I2009bl ) should be multiplied by 2.62 to derive the total surface number density. 



Park et al. 



We display the total surface number density profiles, N C \(R), for the combined, blue and 
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red GCs in Figure [6j We fit the surface number density profile with the projection of the 
NFW profile, n c \(r) = no(r/6) _1 (l + r/b)~ 2 and with the projection of the Dehnen profile, 
n c\( r ) — n o{ r / a ) ~ 7 (1 +r/a) 7-4 . The surface number density profile, N C \(R), is derived from 
the integration of the three dimensional density profile n c \(r) as follows: 



N d (R)=2 r ncl (r)^i= . (3) 

JR 



Vr 2 - R 2 



The solid and long dashed lines represent the projected best fit curves of the NFW 
profile and of the Dehnen profile, respectively. The fitting results for the combined (C) GCs, 
blue (B) GCs, and red (R) GCs are summarized as follows: 

n%(r) = 0.57kp C - 3 (r/6.17kpc)- 1 (l + r/6.17kpc)- 2 

n B( r ) = O.OTkpc-^r/Kmkpc)- 1 ^ + r/10.71 kpcT 2 (4) 

n R( r ) = 0.71kpc~ 3 (r/4.39kpc)- 1 (l + r/4.39kpc)~ 2 , 

for the NFW profile, and 

ng(r) = 0.14kpc" 3 (r/15.87kpc)- 108 (l +r/15.87kpc)- 2 - 91 

n%(r) = 0.03kpc _3 (r/21.14kpc)-°- 94 (l +r/21.14kpc)- 3 - 06 (5) 
n ^f r "^ — n A*7ir-^^-3(^ /q Arl -\^„\-o.70(^ i ^ /q /in 1^^.^—3.30 

for the Dehnen profile. 



0.47kpc" 3 (r/8.40 kpc)-°- 70 (l + r/8.40 kpc) 



It is found that the scale length b of the red GCs in the NFW profile is more than twice 
smaller than that of the blue GCs. This shows that the red GCs are more concentrated 
toward the galaxy center than the blue GCs. 



3.3.2. An Extended Dark Matter Halo in NGC 4636 

We investigate the existence of an extended dark matter halo comparing the velocity 
dispersion profile expected from the stellar mass profile with the measured profiles of stellar 
and GC velocity dispersions for NGC 4636. We first derive the stellar mass profile from the 
surface brightness profile. Then we obtain the velocity dispersion profile expected from the 
stellar mass profile, which will be compared with the measured dispersion profiles for the 
GCs. 



In the left panel of Figure [?l w e plot the surface b rightness profile of NG C 4636 derived 
from o ur KPNO TV - band images (IPark et al.l l2009bl ) compared to those in iPeletier et al. 
(119901 ): iDirsch et al.l (120051 ) for the i?-band photometry. We converted T\ photometry of 
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Park et ah] (l2009bl ) to Cousins i?-band photo metry using th e relati on giv en bvlGeisle r (1996). 
It is seen that the profiles from this study, IPeletier et al.1 (1990|), and iDirsch et al.l (120051 ) 
agree well over the radius. 



We fit the surface brightness profile derived from the KPNO images ( Park et al. 



with the projection of three dimensional luminosity density profile used in 
which is represented by, 



Cote et al. 



2009b) 



(120031 ) 





(T \ V 2 


r 









2( 7 -4) 



(6) 



We obtain, from the fit, the values for the parameters, 7 = 0.84, L tot = 5.42 x 10 10 Lr )Q , 
and a = 2.61 kpc, and overlay the projected best fit curve in Figure [3 We also derive a 
value for the effective radius of R e s = 2/58 ~ 10.99 kpc, which is lar ger than th a t from a fit 
(Reft = 1.'49 ~ 6.347 kpc at Tj-band) using a de Vaucouleurs law in iPark et al.l (l2009bl ). 



In the right panel of Figure [7J we show a three dimensional stellar mass density profile, 
p s (r) = T j(r), with i?-band mass-to-light ratio T = 9.0 M & L~^ & (to be discussed later 
in this section). From the three dimensional stellar mass density profile we derive a stellar 
mass profile of NGC 4636 that is represented by 



MJr) 



Attx 2 p s (r)dx = T / Anx 2 j(x)dx 



'0 

TnZ 



0-^tot 



" (r/a) 1 ' 2 ' 


2(3-7) 


"(7-2 7 ) + (r/a) 1 /2 


l + (r/a)V 2 _ 




1 + (r/a) 1 / 2 



(7) 



We used this stellar mass profile to determine the velocity anisotropy for the NGC 
4636 stellar system and to test the existence of an extended dark matter halo. If we take 
M tot (r) = M s (r) and substitute n c \(r) by p s (r) oc j(r), then we can compute the radial 
component of velocity dispersion profiles (VDPs) of the stars through the Jeans equation 
by assuming the values for _R-band mass-to-light ratios (To) and velocity anisotropics of the 
stellar system [/3 S (?")]• We therefore obtain the projected VDPs for the stellar system from 
the radial component of VDPs through the equation (11) in Section 3.3.4. 

In Figure [8 , we pl ot the stellar VDP for the inner region (R < 34") of NGC 4636 given 



111 



Bender et al.l (119941 ) and the GC VDP derived in this study. The stellar VDP is almost 
constant around 200 km s _1 in the inner region, and is smoothly connected to the GC VDP 
at R ~ 2 kpc that is increasing as the radius increases. Note that the GC VDP varies much 
more than the stellar VDP. In Figure El we also show the projected VDPs calculated with 
T = 8.0 M Q L R l Q} (3 s (r) = 0.0 (isotropic), and T = 9.0 M L^ , (3 s (r) = -0.3 (tangentially 
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biased), which are the best fit curves for the stellar kinematic data of iBender et al.l (119941 ) 
at R < 1 kpc. These values for the mass-to-light ratio and velocity anisotropy for the stellar 
system are si milarly found in the "be st-fitting" halo model using I?-band photometry of 
NGC 4636 bv iKronawitter et all (boooh . 



For the comparison, we also present the projected VDPs calculated using the same 
stellar mass profile as above, but for the GC number density profile n c \(r) and for /3 c i(r) = 
+0.99 (radially biased; upper long dashed line), -99 (tangentially biased; lower long dashed 
line), 0.0 (isotropic; short dashed line). Interestingly, none of these models can account 
for the observed VDPs for the GCs at R > 3 kpc, indicating that mass-to-light ratio is not 
constant over the galactocentric distance, but should increase as the distance increases. This 
demonstrates that there exists an extended dark mater halo in the outer region of NGC 4636. 



3.3.3. X-ray Mass Profiles 



There are several studies that presented the X-ray mass profiles of NGC 4636 (e.g., 



2003 



Mushotzkv et al.lll994l;lBrighenti fc Mathewdll997l ; lMatsushita et alll 998; Lo ewenstein fc Mushotzky 



Johnson et al. 1 120091 ) , as displayed in Figure [9j 



Brighenti fc Mathews! (1997 ) d erived 



a mass profile of NGC 4636 using Einstein HRI observational data of iTrinchieri et al.l (119861 ) 
and ROSAT PSPC observational data of ITrinchieri et al.l (119941 ) . They presented a func- 
tional form of the mass profile, which is used for the an alysis of velocity anisotr opy in this 
study. Their mass profile agrees well with that given by iMushotzky et al.l (119941 ) who used 
ASCA SIS data. 



Loewenstein fc Mushotzky I (120021 . 120031 ) determined a mass profile of NGC 4636 using 
Chandra and XMM-Newton X-ray data. They showed that the total mass of NGC 4636 
increases as r 12 to a good approximation from 0.7 to 35 kpc with ~ 1.5 x 1O 12 M at the 



outermost point. Recently I Johnson et al. I (120091 ) derived a mass profile fo r the inner region 



of NGC 4636 using the Chandra data. It agrees well with other pro files (IMushotzky et al. 



19941 ; iBrighenti fc Mathews! Il997t lLoewenstein fc Mushotzkv 1 120031 ) for 4 < r < 40 kpc, 
while it is much steeper than the others for r < 4 kpc. 



Matsushita et al.l (119981 ) also derived a total mass profile of NGC 4636 using ASCA GIS 



observations, but they did not present the information on the parameters of the mass profile. 
So we read the mass profile in their Figur e 3 and re-derive the mass profile in the functional 
form following the approach of ICote et al.l (120031 ) for M49. ICote et al.l (120031 ) considered that 
the total mass M tot (r) interior to any radius is the sum of dark matter mass Md m (r) (eq. ED 
and stellar mass M s (r) (eq. |7|). If the dark matter density profile is the NFW halo model, 
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Pdm(^) = K(r/r s ) x (l + r/r s ) 2 , the mass profile for dark matter halo is represented by 



M dm (r) 



Atxx 2 pa m (x)dx 



'0 



In ( 1 + - 



(r/r s ) 



1 + (r/r s )_ 

where is the dark matter density normalization and r s is a scale length. 

Fixing M R (r) wi t h ma ss-to-light ratio T = 9.0 M & L~j^ &1 we fit the mass profile given 
by Matsushita et al.l (119981 ) with the total mass profile that is the sum of dark matter mass 
Md m (r) (eq. |SJ) and stellar mass M s (r) (eq. [7]). Thus we determine K = 4.01 x 10 5 M Q kpc -3 
and r s = 147 kpc. 

Then the final mass models of NGC 4636 are 



M tot (r) 
M s (r) 

Mdm(r) 



MJr) + M dm (r) 



4.87 x 10 U M P 



(r/2.61kpc 



,V2 



1 + (r/2.61 kpc) 1 / 2 
5.33 + (r/2.61 kpc) 1 / 2 



4.33 



X 



1 + (r/2.61 kpc) 1 / 2 
il x 1O 13 M 

In I 1 + 



(r/147kpc) 
147kpc7 (l+r/147kpc) 



(9) 



In Figure [H we plot the stellar mass profile (dotted line), the dark matter halo (long- 
dashed line), and the total mass profile (dot-dashed line) derived from this method. 



While four mass profiles given bvlMushotzkv et al.l (119941) ; iBrighenti fc Mathewd (119971 ) ; 
Loewenstein fc Mushotzky I (120031 ) ; I Johnson et ah 1 J2009h show monotonic increase with in- 
creasing radius, only the mass profile given by iMatsushita et al.l (119981 ) shows a flattening 
at r m 10 — 20 kpc, and then keeps increasing thereafter out to r « 400 kpc. From this 
compa rison we conclude that the mass profile for r > 10 kpc given by IMatsushita et al. 
( 119981 ) may be in error. Since X-ray mass profiles that account for the dark matter halo are 
similar to or larger than the stellar mass profile at r > 2 kpc, we consider that the X-ray 
mass profiles are good enough to determine the velocity anisotropy of the GC system for the 
following analysis. 
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3.3.4- Determination of the Velocity Anisotropy 



The velocity anisotropy of GCs is determined as follows. First, assuming the velocity 
anisotropy [/? c i(r)] in prior, we derive the theoretical projected VDP [a p (R)] and theoretical 
projected aperture VDP [a ap (< R)] using the Jeans equation. The projected aperture VDP is 
the velocity dispersion of all objects interior to a given projected radial distance R. To obtain 
those theoretical VDPs, we use the GC number density profile [n c i(r)] of the combined, blue, 
and red GCs and the mass profile [M tot (r)] derived in the previous section. Second, from 
the comparison of these calculated VDPs with measured VDPs, we determine the velocity 
anisotropy of GCs. 

We begin by deriving the theoretical projected VDPs. The spherical Jeans equation 
(eq. [2]) can be solved for the radial component of velocity dispersion, o~ r (r): 



a 2 (r) 



n c i(r) 



x 



exp 



2& 



dr 



GM tot 
n d — exp 

x A 



2/3, 



cl 



X 



dx dx 



(10) 



Then the projected VDP, cr p (R) can be derived by 



N cl (R) 



n cl a 2 r (r) 1-/3, 



R 2 



r dr 



Vr 



R 2 



(11) 



The projected aperture VDP, cr ap (< R), can be computed by 

r r R 

o-l p (<R) = / N cl (R')a 2 p (R') R' dR' 

_»/i? m i n 

x / N cX (R')R'dR' 

^min 



(12) 



where -R m i n is the projected galactocentric distance of the innermost data point in the GC 
sample (-R m i n = 1.63 kpc in this study). 

We present the measured VDP in comparison with the VDPs calculated by assuming 
several velocity anisotropies in Figures [TUl and [TT1 The left and right panels in Figure [TU1 show 
the VDPs calculated using the NFW profile and the Dehnen for the GC number density, 
respectively. The upper panels show the projected VDPs, and the lower panels show the 
projected aperture VDPs. The measured dispersion data taken from Figure [5] are shown by 
filled circles along with their confidence intervals. The projected aperture VDPs in the lower 
panels are plotted in the similar fashion to the case of the upper panels. 
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The mean velocity dispersio n for the outer par t show s a difference by a factor of ~ 1.5 be- 
tween the VDP profiles based on IMatsushita et al.l (119981 ) and those on lBrighenti &: Mathews 



( 119971 ) and lLoewenstein &: Mushotzkv I (120031) . This is because the mass in the o uter part 



derived from Brighenti fc Mathews! (19971) an d Loewenstein fc Mushotzkv ( 2003 ) is about 



twice larger than that from iMatsushita et al. ( 19981 ). As a result, the sys tem of all the GCs 
has a t ange ntially biased velocity ellipsoi d (0 r ] < 0) for the mass profiles by lBrighenti fc Mathews 
( 119971 ) and lLoewenstein fc Mushotzkv I (120031 ). but has a radially biased velocity ellipsoid 
((3 c i > 0) for the IMatsushita et al.l (119981 ) mass profile, in Figure (TU] (a,b) based on the NFW 
profile for the GC number density. A similar result can be found in Figure [10] (c,d) based 
on the Dehnen profile for the GC number density. We tried an eye-ball fit to the data for 
the outer region of NGC 4636, find ing that the VDP profile with (3 C \ = —9.0 (derived for the 
Loewenstein fc Mushotzkv I (120031 ) mass profile) fits approximately the data for R > 20 kpc 
(shown by the solid lines labeled with f3 c \ = —9.0 in (b) and (d)). 

In Figure [TTJ we show a similar analysis for the blue (left panels) and red GCs (right 
panels) based on the NFW profile for the GC number density. These figures show similar 
tre nds to those for all the GCs. T hey have a tangentially biased ve l ocity ellipsoid (j3 c \ < 0) for 



the lBrighenti fc Mathewsl (119971 ) and lLoewenstein fc Mushotz! 



kv 



(boo3h mass p r ofiles, while 



they have a radially biased velocity ellipsoid {(3 C \ > 0) for the IMatsushita et al.l (119981 ) mass 
profile. When we compare the velocity anisotropies between the blue and red GCs by fixing 
the m ass p rofile, we find the velocity anisotrop ies of the red GCs by the lBrighenti &: Mathews 



( 119971 ) and lLoewenstein fc Mushotzkv 1 120031 ) mass profiles are slightly more tangential than 
those of the blue GCs in Figure [TTJ These trends are similar to the velocity anisotropies 
based on the Dehnen profile for the GC number density, too, although they were not plotted 
in Figure [TTJ We tried an eye-ball fit to the dat a for R > 20 kpc, finding that the VDP 
profiles with (3 C \ = —1.0 and -25 (derived for the lLoewenstein fc Mushotzkv I (120031 ) mass 
profile), fit approximately the data for the blue GCs and red GCs, respectively (shown by 
the solid lines labeled with (3 C \ = —1.0 and -25, respectively, in (b) and (d)). It is noted 
that the velocity dispersion for the blue GCs increases as R decreases so that the orbit of 
the blue GCs becomes radial in the inner region at R < 12 kpc. 



Since we concluded that the mass profile given by IMatsushita et al. 



ror be fore, we adopt finally the results based on the mass profile giv e n by 



( 199 81) may be in er- 



Loewenstein fc Mushotzkv 



( 120031 ) (that is also similar to that given bv lBrighenti k Mathewsl jl997h ). [n summary, it is 



found that the orbit of the GC system in NGC 4636 is tangential. The orbits of both the red 
GCs and blue GCs are tangential, and the orbits of the red GCs is slightly more tangential 
than that of the blue GCs. 
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Discussion 



4.1. GC Kinematics and the Global Properties of gEs 

In this section, we combine the kinematic properties of the GC systems in gEs by 
including the results for the N GC 4636 GC system i n thi s study, those for the NGC 1407 
GC system (172 GCs) given i n IRomanowskv et al.l (120091). and those for the GC systems 



in other gEs in iHwang et al. 
276 G Cs in M87 JCote et all 



( jLee et al 



1399 faichtler et al.l l2004h7341 GCs in NGC 512 8 JWoodlev et 



2008a|), and 172 GCs in N GC 4636 flSchuberth et al 



method adopted in this study (see also IRomanowskv et al 



(2008|). IHwang et al.l (120081) analyzed t he velocity data of 



200 lh . 263 GCs in M 49 (lC6ted~al1 I2QO3J) 435 GCs in NGC 

J2003), 121 GCs in M60 
20061 ). using the similar 



20091 ) . We de rived the kinematic 



param eters of the NGC 1407 GCs from the data in the catalog given by IRomanowskv et al. 



( 120091 ) using the same analysis as used for other gEs in this study, and listed the results in 
Table [2j The information on the host galaxies is listed in Table [3X and the global kinematics 
of GCs in gEs are summarized in Table HI 



The velocity dispersion of the GCs in NGC 4636 is similar to that of the M60 GCs 
[a p = 234^ km s _1 ), although NGC 4636 is about one magnitude fainter than M60. This 



indicates that the mass to luminosity ratio is larger in NGC 4636 than in M60. The velocity 
dispersion of the GCs in NGC 4636 is much smaller than those for the three brightest gEs 
(M87, M49 and NGC 1399). However, the velocity dispersion of the GCs in NGC 4636 is 
about twice larger than than that of NGC 5128 GCs (a p = 129iy km s" 1 ), although NGC 
4636 is slightly fainter than NGC 5128. 

The rotation-corrected velocity dispersion, a v , r of the blue GCs in NGC 4636 is larger 
than that of the red GCs. This trend is similarly seen in the three X-ray brightest gEs. 
However it is opposite in the case of M60, and the rotation-corrected velocity dispersions 
of both blue and red GCs are similar in NGC 5128 and NGC 1407. The rotation-corrected 
velocity dispersion ranges from 129 km s" 1 (NGC 5128) to 399 km s" 1 (M87) for all the 
GCs, and similarly for the blue GCs and red GCs. The rotation amplitude ranges from 30 
km s' 1 (NGC 5128) to 172 km s" 1 (M87) for all the GCs, and similarly for the blue GCs and 
red GCs. The ratio of the rotation amplitude to the rotation-corrected velocity dispersion 
ranges from 0.10 (NGC 1399) to 0.65 (M60). M87 and M60 show much stronger rotation 
with VLR/GpG c > 0.4 than the others. 

We investigate any dependence of the kinematic properties of the GC systems in gEs 
on the global properties of their host galaxies. We calculated the Spearman's correlation 
coefficient (r^) and its significa nce (a(rs)) to ch eck any correlation between parameters, 
and used the bisector method (llsobe et al.l Il990l ) to do linear fits for correlated pairs of 
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parameters. 

Figure [redisplays three kinematic parameters (the rotation-corrected velocity dispersion 
for all the GCs, &p^ c (as well as the blue and red GCs), the ratio of the rotation-corrected 
velocity dispersion between the blue GCs and red GCs, &pfi° / &p^° , and the ratio of the 
rotational velocity to the rotation-corrected velocity dispersion for all the GCs, flR/dp^ ') 
versus three global parameters (X-ray luminosity, log L x , stellar velocity dispersion, a star , 
and total l^-band magnitude, M v ) for the gEs, where the rotation-corrected velocity dis- 
persion means the velocity dispersion about the best fit rotation curve. Here the stellar 
velocity dispersion represents the mean value o f the velocity stellar dispersion at ~ -R e ff/4 



for gi ven galaxy: NGC 1 399 flSariia et al.ll2000h. NGC 5128 flWilkinson et al.lll986f ). NGC 



1407 (jSpolaor et al.1 120081 ) and other galaxies (IBender et all 11994 ) 



The rotation-corrected velocity dispersion for all the GCs shows a strong correlation 
with the X-ray luminosity and stellar velocity dispersion, and a weaker correlation with the 
total V-band magnitude, as seen in Figure H2T a). (b) and (c). This shows that the rotation- 
corrected velocity dispersion for the GCs is an excellent indicator for the mass or luminosity 
of their host galaxies. It is noted that the stellar velocity dispersion of M60 is as large as 
that of M87 (see Fig. [T2T b)). being much larger than that expected from its luminosity. This 
may be related with the presence of companion SBc galaxy NGC 4647 located at 2''. 5 from 
M60, but the cause of this large stellar velocity dispersion of M60 is not known. M60 seems 
to be different from the other gEs in stellar velocity dispersion. Therefore we excluded any 
parameters for M60 related with the stellar velocity dispersion for the following correlation 
analysis. 

Linear fits to the data yield: = 97.66 log Lx — 3796 with rms = 28 km s _1 , 

cjAgc = i.394o- stor _ 69.51 with rms = 31 km s~ l , and a^ c = -202.8M V - 4201 with rms 
= 76 km s _1 . In the case of a Ptr and logLx, the blue GCs show a stronger correlation than 
the red GCs: <Jp^ c = 110.28 logLx — 4309 with rms = 32 km s _1 and the Spearman's rank 
correlation coefficient rs = 0.96 for the blue GCs, and &p^ c = 84.04 log — 3243 with 
rms = 27 km s _1 and rs = 0.86 for the red GCs. It is noted that the slope for the blue 
GCs (110.28) is steeper than that for the red GCs (84.04). On the other hand, it is opposite 
in the case of a p>r and a star . The red GCs show a stronger correlation than the blue GCs: 
a^ c = 1.237<w - 50.31 with rms = 24 km s" 1 and r s = 0.0.83 for the red GCs, and 
< 7 p,r !C = l-538cr star — 88.77 with rms = 46 km s _1 and rs = 0.71 for the blue GCs. It is noted 
that the slope for the blue GCs (1.538) is steeper than that for the red GCs (1.237). 

The slope for the relation between (Jp^ c and a star (solid line) is larger than one (dotted 
line). This indicates that the more massive gEs are, the more massive dark matter halo 
they have. The relation between cr^f and My shows a larger scatter than that between 
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(Tp^ c and X-ray luminosity. However, NGC 1399 has an absolute magnitude much fainter 
than expected from its X-ray luminosity. If NGC 1399 is removed, the linear fit yields 
°pr = ~ 255. 6My — 5389 with a smaller scatter, rms = 65 km s _1 (dashed line in Fig. 
IT27c)). On the other hand, little correlation is seen between other pairs of parameters. 

In Figure [131 we display the same three kinematic parameters versus three other global 
parameters (specific frequency, Sn, the number ratio of the blue GCs and the red GCs, 
Nbgc/Nrgci an d ellipticity, e) for the gEs. We find strong correlations for some pairs 
of parameters: (a) a^ c and S N , (b) a^ c and e, (c) a^GC ' / a RGC and N BGC /N RGC , and 
(d) VtR/a^ c and N BG c/N RGC . Linear fits to the data yield: a^ c = 20.03SW + 150 
with rms = 63 km s _1 (panel a), &pf c = — 1548e + 552 with rms = 57 km s" 1 (panel 
c )> <7p,? C /<7p? C = -0.90SN BGC I 'N RGC + 2.16 with rms = 0.11 (panel e), and VLR/a^ c = 
1.120N BGC /N RGC -1.00 with rms = 0.09 (panel h). NGC 1407 shows much different relations 
including N BGG /N RGG from those of other gEs so that we did not use NGC 1407 for linear 
fitting for (c), (e) , and (h). Thus the rotation-corrected velocity dispersion for all the 
GCs increases as Sn increases, but this velocity dispersion decreases as ellipticity increases. 
a p,r !C / a p^ C decreases as N BGG / N RGG increases, while (flR) / (T^ c increases as N BGG /N RGG 
increases. 

Figure [T4T a) displays directly the relation between the velocity dispersion and rotational 
amplitudes of the GCs: The ratio of the velocity dispersion between the blue GCs and red 
GCs [pp^ c l&p^ c ) vs. the ratio of the rotational velocity to the velocity dispersion for all 
the GCs {Q.R/a^ c ). It is found that a^ c /cr^ c has a strong correlation with VtR/a^ : 
VtR/a^ c = -\.\$&t*° c I<t™* c + 1.556 with rms=0.108. That is, the weaker the rotation 
is, the larger the velocity dispersion ratio between the blue GCs and red GCs is. 

Figure [T4Tb) displays directly the relation between the rotational amplitudes and rota- 
tional axis of the GCs: The ratio of the rotational velocity to the velocity dispersion for all 
the GCs (flR/dp^ ) vs. the difference between the GC rotation angle and the position angle 
of the minor axis of their host galaxies (Bo — Qminor)- Interestingly it is seen that the GC 
rotation axis for the two gEs with the strongest rotation (M60 and M87) is approximately 
aligned with the minor axis of their host galaxies (the dotted lines in the shaded regions). 
This shows that the GC systems in gEs rotate around the minor axis of their host galaxies, 
once they have some strong rotation. 

In Table [5j we summarized the results of the linear fits (Y = aX + b) given in Figures 
[T2l [TBI and [TH We listed the Spearman's rank correlation coefficient (r§) and its signifi- 
cance (a(rs)) in the 6th and 7th columns , respectively. These correlations found in this 
study provide strong constraints on modeling the GC systems in gEs (to be discussed in 
§4.3p . Moreover, it is needed, in the future studies, to investigate which relation is the most 
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fundamental one, if any, among several correlations between kinematic properties of the GC 
systems in gEs and the global parameters of their host galaxies found in this study. 

Figure [TBI represents the velocity anisotropy versus the velocity dispersions for the GCs 
in 6 gEs (M6 M87 , M4 9 NGC 1399, NGC 1407, a nd NGC 4636), determined in this study, 
Hwang et all (120081 ) and iRomanowskv et al.1 (120091 ). We consider [3 C \ = 0.5, 0, -1, and -99, 
as the radial, isotropic, tangential, and strongly tangenti al velocity anisotropy, respe c tively . 
In the case of NGC 4636, we used the result derived for lLoewenstein fc Mushotzky I (120031 ) 
mass profiles. All the GC systems in three gEs (M49, M87 and NGC 1399) show isotropic 
orbits, while those in two gEs (M60 and NGC 4636) show tangential orbits. It it noted also 
that more massive gEs (with larger GC velocity dispersion) have isotropic orbits, while less 
massive gEs have tangential orbits. The blue and red GCs in three gEs (M49, NGC 1399 
and NGC 4636) show similar orbits, while those in two gEs (M60 and M87) show opposite 
orbits. 



4.2. Radial Variation of the GC Kinematics for gEs 

To investigate the radial variation of the kinematics of the GC systems in the gEs, we plot 
the rotation-corrected velocity dispersions {o- Jhr / a^ c ) against the projected galactocentric 
distances in Figure [HJ The rotation-corrected velocity dispersion is normalized by that of 
all the GCs in each gE. The projected galactocentric distance is normalized with respect 
to the effective radius of each gE, R/R c g. The combined GCs do not show any significant 
change of the velocity dispersion over the whole region of a galaxy. However, the mean 
velocity dispersion for the red GCs in the inner region (R < 2R e g), 1.05 ± 0.15, is slightly 
larger than that in the outer region (R > 2R e g), 0.85 ± 0.07, while that for the blue GCs in 
the inner region is little different from that in the outer region. Note that the red GCs in 
M60, NGC 5128, and NGC 1407 have much larger a v , T jo^ c in the inner region than in the 
outer region. Why the red GCs have larger velocity dispersion in the inner region than in 
the outer region is an interesting question to solve. 

Figure [T7] displays the ratio of the rotation amplitude to the rotation-corrected velocity 
dispersion {VtR/ a^ c ) vs. the projected galactocentric distance normalized to the effective 
radius (R/R cS ). The combined GCs in M87, NGC 5128, and NGC 1407 show much larger 
VtR/a^ c in the outer region than in the inner region, while those in other galaxies show 
little, if any, radial variation. The blue GCs and red GCs in M87 and NGC 1407 show a 
similar trend, while only the red GCs in NGC 5128 show a similar trend. However, in the 
case of NGC 4636, both the blue and red GCs show an opposite trend, a larger QR/a^ 
in the inner region than in the outer region. It is noted that the blue GCs and red GCs in 
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NGC 4636 rotate in opposite direction to each other. 

To show better the results seen in Figures [H] and [T71 we display in Figure [TS] the ratio of 
the rotation-corrected velocity dispersion between the blue GCs and red GCs (&p ^ c '/ >cr p^ C ) 
and the ratio of the rotation amplitude to the rotation-corrected velocity dispersion between 
the blue GCs and the red GCs ( (QR) / ap 3 ^ / (QR) /a^ ) as a function of projected galac- 
tocentric distance normalized to the effective radius (R/Rgg). It is seen that the mean value 
of &p t r C / a pf C f° r the inner region (R < 2_R eff ), is close to one, 0.97 ± 0.19, and that it is 
slightly smaller than that for the outer region (R > 2i? eff ), 1.26 ± 0.12. It is also noted that 
the values of (QR) / <Jpf c / (QR) / u^ c for the inner region, 1.77 ± 1.34 are marginally larger 
than those for the outer region, 1.14 ± 0.87. These results show that, when compared with 
the red GCs, the blue GCs have a larger velocity dispersion and smaller rotation amplitude 
in the outer region. 



4.3. Comparison with Formation Models of Globular Clusters in gEs 



The kinematic properties of the GC systems in galaxies can pro yide strong constraints 



on the formation models of the GC system s and their host galaxies (ICote et al 



2001, 2003 



Hwang et al.ll2008l ; iKormendy et al. 1120091 ). Several formati on models of the GC systems in 
gEs have been suggested t o describe the formation of GCs ([Peebles I Il969l ; lAshman &: Zepi 



1992 



Harris et al. 1995 ; Forbes et al. 1997 ; Cote et al. 1998h and a summary of model 



descriptions and predictions can be found in several literature (IK 
Richtler et all booi IWest et all booi iBrodie fc Straderi booa 



lode fc Z epf 2001 



Lee 2003 



Hwang et all 120081 ). These 



models can be broadly divided into four categories: the monolithic collapse model, the major 
merger model, the multiphase dissipational collapse model, and the dissipationless accretion 
model. These models suggest some kinematic properties of the GCs in galaxie s. These were 



discus sed in comparison with the observation results for the GCs in gEs in iHwang et al. 



(b()08h and will be summarized below. We focus on kinematic aspects of these models to 
compare with observational results below. It is noted that these classical models do not 
provide detailed kinematic properties so as to compare with the observational results. 

In the monolithic collapse model an elliptical galaxy and its GCs are fo rmed through 
the collapse of an isolated massive gas cloud or protogalaxy at high redshift (jLarsonlll975l ; 
Carlber Jl98i ; lArimoto fc Yoshdll987h. In this m odel, the rotation of GCs can be generated 
by tidal torques from companions (IPeebles Ill969l ). but the resulting rotation is not expected 
to be strong. This model cannot explain the strong rotation of the GC system seen in M60 
and M87, and the globally isotropic velocity ellipsoid of the GC system in several gEs. 



-20- 



In the major mer ger model elliptical galaxies are formed by a merger of two or more 
gas-rich disk galaxies (lToomrelll977l ; lAshman &: ZepJll992l ; IZepf et al.ll2000l ). In this model, 
younger, spatially concentrated, red GCs are formed during the mer ger, while spatiall y more 



extended, blue GCs come from the halos of the disk galaxies (e.g.. lBekki et al.ll2002l ). This 
model predicts that the newly formed red GCs show little rotation compared to the blue GCs 
since the angular momentum would be transported to the outer region during the merging 
process. This model cannot explain the presence of rotation of the red GCs in M60, M87, 
and NGC 5128, while its prediction is consistent with the absence of rotation of the red GCs 
in other gEs. 



In the multiphase dissipational collapse model (IForbes et al.l Il997l ) elliptical galaxies 
form their GCs in distinct star formation phases through a dissipational collapse, and they 
capture some GCs by tidal effects from neighboring galaxies or the accretion of dwarf galaxies. 
The blue GCs are formed in the first star formation phase and the red GCs are formed in 
the subsequent star formation phase after the gas in the galaxy is self-enriched. This model 
predicts that the blue GC system shows no rotation and a high velocity dispersion, while the 
red GC system shows some rotation depending on the degree of dissipation in the collapse. 
Their prediction for the rotation is not consistent with the observational results for gEs in 
this study (e.g., rotation measured in M60 and M87). 



In the dissipationless accretion model (jCote et al.lll998l ) the red GCs are formed in a 
dissipational monolithic collapse of a primary proto-galactic cloud, while the blue GCs are 
subsequently captured from other (low-mass) galaxies through mergers or tidal stripping. 
Since the blue GCs are captured from other galaxies, they are expected to show a spatial 
distribution more extended than that of the red GCs. The blue GCs are expected to have 
radially biased orbits r ather than isotropic or tangentially biased orbits, and are also expected 
to show no rotation (iRichtler et al.ll2004l ). These predictions are not consistent with the 
observational results that the blue GCs in gEs show isotropic or tangential orbits as seen in 
Fig- [15T b). and that the blue GCs in four gEs show measurable rotation as seen in Table HI 

Above classical models give some qualitative predictions for the kinematic properties of 
the GC systems, but without any quantitative information. On the other hand, numerical 
simulations provide quantitative predictions for the kinematic properties of the GC systems 
in galaxies, which can be compared with the observational results. There are several nu- 
merical simulation studies th at provided some pr e dictions of kinematic properties of the 



GC systems in galaxies ( e .g.. IVesperini et al. 



2003 



Bekki et a 



Kravtsov fc Gnedinl [2003 : iMoore et alihood : iBekki et al.ll2008l ). Focusing on kinematic as 



2005; Diemand et al. 2005 



pects of the GC systems in gEs, we compare our observational results for several gEs with 
the results in the simulations. 



-21 - 



Bekki et al.l (120051 ) presented the results of a simulation of dissipationless major mergers 
of spiral galaxies derived with an assumption that the spatial distribution of the GCs in 
E/SO's are initially similar to that for the Milky Way Galaxy. They predicted that both pre- 
existing metal-poor globular clusters (MPGCs) and metal-rich globular clusters (MRGCs) 
obtain stronger rotation in the outer region regardless of the orbital configuration of the 
merging galaxies. In Figure [TTJ, the blue GCs and red GCs in M87 and NGC 1407 show 
slightly stronger rotation in the outer region (R > 2_R eff ) than in the inner region (R < 2i? eff ), 
which is consistent with their prediction. However, the blue GCs and red GCs in other gEs 
are not consistent with their prediction. 



Bekki et al.l (120051 ) also predicted that the velocity dispersion for both the MPGCs and 
MRGCs decreases as the galactocentric distance increases and that there is little difference 
between the MPGCs and MRGCs in all major merger models, while those are sometimes 
flat in multiple merger models. In addition the MPGCs should show slightly larger central 
velocity dispersion than the MRGCs, indicating that the MPGCs are dynamically hotter 
than the MRGCs. However, the observational results for NGC 4636 in Figure 5 and other 
gEs in Figure US show that there is significant difference in the VDPs between the MPGCs 
and MRGCs, which is not consistent with the prediction by iBekki et al.l (120051 ). 



Recently iBekki et al.l (120081 ) investigated the origin of GC systems using high- resolution 
cosmological N-body simulations combined with the semi-analytic models of galaxy forma- 
tion. They predicted that the majority (about 90 per cent) of GCs seen in the galaxy halos 
today were formed in low-mass dwarf galaxies at redshifts larger than 3 and that the mean 
formation epochs for the MPGCs ([Fe/H]< -1) and MRGCs ([Fe/H]> -1) are, respectively, 
12.7 Gyr ago (redshift z = 5.7) and 12.3 Gyr ago (redshift z = 4.7). MPGCs and MRGCs 
were also formed in the gas-rich major mergers for a wide range of redshifts, with peaks at 
z ~ 5 and z ~ 3, respectively. MRGCs are formed also in isolated gas-rich galaxies. They 
could explain several observational results on structural, kinematic and chemical properties 
of the GC systems in various kinds of galaxies with their models. 



Kinematic properties of the GC systems predicted in the IBekki et al.l (120081 ) 's model 
are: (a) the velocity dispersions of the MPGCs and MRGCs increase according to the total 
luminosity of their host galaxy (Mb)', (b) the ratio of the velocity dispersions between the 
MPGCs and MRGCs is almost one for a wide range of the total luminosity of their host 
galaxy; and (c) the ratio of maximum rotational velocities to central velocity dispersions in 
GC systems is low (V/a < 0.3) for most galaxies, which is i n contrast to the c ase of disk-disk 
major mergers leading to larger V/cr (> 0.5) suggested by IBekki et al.l (120051 ). 



Their prediction (a) is roughly consistent with the observational results shown in Figure 
IT2Tc). Their prediction (b) also appears to be consistent with our observational results 
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consi dering the large uncertainties (see Fig. IT27 f) in this study and Fig. 18 in IBekki et al. 
20081 ). However, it is worth noting that the ratio of the velocity dispersions between the blue 
and red GCs ranges from 0.85 to 1.25 in Figure H2T f). that the ratio is the largest for the 
sample of the brightest and faintest galaxies, and that it is the smallest for the intermediately 
luminous galaxy. This trend should be checked with a larger sample of galaxies. Four out 
of seven gEs (NGC 4636, NGC 1399, M49, and NGC 5128) have low ratios of maximum 
rotational velocities to central velocity dispersions (V/cr < 0.3), as shown in Figure fTWg). 
which is again in a broad agreement with their prediction (c). 

Therefore, the o bservational r esults appear to be approximately consistent with those 
in the simulations of IBekki et al.l (120081 ) in terms of the parameters that they considered. 
However, it is noted that some GC systems show different behavior depending on the pa- 
rameters and there remain many observational results to be explained (e.g., Figs. [T21 and 

nm 



4.4. A Mixture Scenario for the Origin of Globular Clusters in gEs 



The most notable results emerging from our study of kinematic properties of the GC 
systems in seven gEs are (a) that the kinematic properties of the GC systems are diverse 
among gEs and (b) that some kinematic parameters of the GC systems show strong cor- 
relations with the global parameters of their host galaxies. The first result indicates that 
the GCs in gEs were probably formed and evolved via diverse mechanisms rather than via 
one single way, and the second result implies that the kinematics of the GCs is controlled 
by the gravitational potential of their ho st galaxie s . Considering all above comparisons and 
observational aspects of the GCs in gEs (lLeell2003c iHwang et al.l 12008c IPeng et al.l 120081 ). we 
derive a following scenario for the origin of GCs in gEs. 

(1) MPGCs are formed mostly in low-mass dwarf galaxies very early, and preferentially 
in dwarf galaxies located in the high density environment like galaxy clusters. These are the 
first generation of GCs in the universe. MPGCs should be also formed in massive galaxies 
as well, but the number of these massive galaxies is much smaller compared with that of 
the dwarf galaxies. Observational evidence for this are: (a) MPGC are among the most 
metal-poor objects; (b) Th e spatial distribution of the MPGCs is more extended t han that 
of the MRGCs around gEs Lee. Kim, fc Geisler E99I birsch etalll2005cEe3l2003h : (c) The 
dwarfs in higher densit y environment ha ve higher specific frequency than those in the lower 
density environments (IPeng et al.l 120081 ); and (d) the intracluster GCs are mostly MPGCs 



flTamura et al. 112009 ) 
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(2) MRGCs are formed together with stars in massive galaxies or dissipational merging 
galaxies later than MPGCs, but not much later than MPGCs. The chemical enrichment of 
the galaxies is rapid after the formation of MPGCs and the difference in the formation epoch 
of MPGCs and MRGCs should be small (e.g., much smaller than 1 Gyr). Observational 
evidence for this are: (a) MRGCs have, on average, sim ilar metallicity to that of the stellar 
halo in gEs (jGeisler. Lee, fc Kim Hl998uLee et al.ll2008bl ); and (b) the difference in estimated 
ages between the MPGCs and MRGCs is s mall, whil e the difference in mean metallicity of 
the two populations of GCs is about a dex (jLeell2003l ; iLee et al.ll2008bl ). 



(3) Massive galaxies grow becoming gEs via dissipationless or dissipational merging of 
galaxies of various types and via accretion of many dwarf galaxies. New MRGCs will be 
formed during dissipational merging, but the faction of dissipational merging at this stage 
should be minor. A significant fraction of MPGCs in gEs we see today are from dissipationless 
merging or accretion. Observational evidence for this are: (a) The s patial distribution of 



the MPGCs is more extended than that of the MR GCs around gEs (ILee. Kim, fc Geisler 



1998; Lee 



2003: 



MPGCs ( ITamura et al. 



Dirsch et al.ll2005l ; ILee et al.ll2008bl ): (b) The intracluster GCs are mostly 
2006); (c) The kinematics of the GCs are diverse among the gEs, 



and there are strong correlations of the GC kinematics with some global p a ramet ers of their 



host galaxies found in this study. It is also noted that iKormendy et al. I (120091 ) concluded 



from the study of structure of a large sample of elliptical galaxies that bright boxy elliptical 
galaxies (like gEs) were formed via dissipationless (dry) merger, while faint disky elliptical 
galaxies were formed via dissipational (wet) merger and conversion of late-type galaxies into 
spheroidals. 

In this scenario each gE has a different history of growing involved with diverse merging 
and accretion, explaining naturally the diversity in the kinematics of the GC systems in 
gEs. This scenario also explains the bimodal color distribution of the GCs, the difference in 
spatial distribution between the blue GCs and red GCs, the correlation in color between the 
red GCs and their host galaxies. This is a mixture model or a bibimbap model (bibimbap is 
a Korean dish where diverse vegetables and warm rice are mixed together with some tasty 
sauces) in that it includes all the key elements in the previous models. 



5. Summary 

Using the photometric and spectroscopic database of 238 GCs (108 blue GCs and 130 
red GCs) in NGC 4636 (gE in the Virgo cluster), we have derived the kinematics of the GC 
system of this galaxy. Then we have compared the kinematics of the GC systems in seven 
gEs including NGC 4636, and have investigated correlations between the GC kinematics 
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and the globular parameters of their host galaxies. Our primary results are summarized as 
follows. 

1. The red GC subsample of NGC 4636 shows marginal overall rotation. In the inner 
region both the blue and red GC subsamples show some rotation. 

2. Both of the velocity dispersions about the mean velocity and about the best fit rotation 
curve of the blue GCs are about 50 km s _1 larger than those of the red GCs. 

3. Comparison of observed stellar and GC velocity dispersion profiles with those calcu- 
lated from the stellar mass profile shows that the mass-to-light ratio is not constant, 
but should increase as the galactocentric distance increases, indicating the existence of 
an extended dark matter halo in NGC 4636. 

4. Using the X-ray mass profiles, the number density distribution of GCs, and the ob- 
served VDP of GCs, we have determined the velocity anisotropy of the NGC 4636 GC 
system. The orbits of the NGC 4636 GC system are tangentially biased. The orbits 
of both the red GCs and blue GCs at R > 20 kpc are tangential, and the orbits of the 
red GCs is slightly more tangential than that of the blue GCs. The orbit of the blue 
GCs in the inner region at R < 12 kpc is found to be radial. 

5. From the comparison of the kinematics of the NGC 4636 GC system in this study with 
those for other gEs, we found several correlations between the kinematic properties of 
the GC systems and the global parameters of their host galaxies. 

6. We compared the observational results for the GC systems in gEs including NGC 4636 
with those in several GC formation models, and found that some results are consistent 
with the predictions from the models but some are not. We suggested a mixture 
scenario for the origin of the GCs in gEs. 
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Table 1. Kinematics of the NGC 4636 Globular Cluster System 
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Table 2. Kinematics of the NGC 1407 Globular Cluster System 
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217 
666 


223 
118 
357 


86 
48 
38 


175011 

17231^ 
17811g 


25llg 
278!*; 
214l| 


2091^ 
990+40 

zzo -50 

18911 


10411 
8011 
15011? 


250l^ 7 
909+39 

1991? 7 


41 +a17 
0.281°;! 
7S+ - 29 

u - ' °-0.25 
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Table 3. Giant Elliptical Galaxy Samples 



Galaxy 


M v a 


v b 

u sys 


Reff C 


e d 


PA ■ e 
l ..n.miii 


Distance * 


O'star 9 


log(ix) h 


N G c 1 


Sat j 






(km s" 1 ) 


(kpc) 




(deg) 


(Mpc) 


(km S _1 ) 


(erg s" 1 ) 


blue red 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


M60 


-22.13 


1056 


7.33 


0.216 


15 


16.8 


337± 9 


41.27±0.042 


663 455 


4.0±0.8 


M87 


-22.38 


1307 


7.16 


0.125 


69 


16.1 


330± 5 


43.08±0.007 


928 771 


14.1±1.5 


M49 


-22.57 


997 


9.48 


0.175 


65 


16.3 


294± 5 


41.71±0.043 


683 563 


3.6±0.6 


NGC 1399 


-21.71 


1442 


14.55 


0.099 


20 


20.0 


250±12 


42.18±0.040 


459 500 


5.1±1.2 


NGC 5128 


-21.90 


541 


6.02 


0.224 


125 


4.2 


140±40 


40.15±0.200 


178 158 


1.8±0.5 


NGC 4636 


-21.67 


928 


6.35 


0.256 


58 


14.7 


202± 3 


41.68±0.046 


639 633 


8.2±1.6 


NGC 1407 


-21.86 


1784 


7.34 


0.050 


148 


21.0 


265±15 


41.14±0.092 


446 725 


3.8±1.3 



a V-band absol ute total magnitude: NGC 5128 jPufour et al.l Il979l l. NGC 1407 l|Bender et all Il992h . and other galaxies 
l lFaber et a,l.lll997tl . 

b S ystemic velocity: M60 llLee et alj|2008blj. M87. M49 llSmith et al lbOOOh . NGC 1399 jflichtler et alj|2004l . NGC 5128 llHui et al.l 
Il995h . NGC 4636 (Paper I), and NGC 1407 jRomanowskv et alJl2009h . 

c Effective r adius in units of kpc: M60 llLee et alj|2008bh . NGC 5128 l lDufour et al.lll979h , NGC 1407 l lSpolaor et al.ll2008Tl , and 
other galaxies l|Kim et al.ll2006T l. 

d Ellipticity: NGC 5128 llde Vaucouleurs et 3.1.^199^ . NGC 1407 l lSpolaor et al.ll2008h . and other galaxies llKim et al.ll2006l) . 

e P osition angle of the minor axis: M60 <Lee et al.ll2008bh. N GC 1399 jSaelia et alj|200ch . NGC 5128 llDufour et al.lll979h . NGC 
1407 llSpolaor et alj|2008h , and other galaxies llKim et al.ll2006l) . 

f Distance in units of Mpc: NGC 1407 l lSpolaor et al.ll2008l) . and other galaxies Flbnrv et all ll200ll l. 

g Mean velocity s tellar dispersion at ~ R eB /4: NGC 13 99 llSaglia et alj|200d) . NGC 5128 llWilkinson et al.lll986h . NGC 1407 
llSpolaor et alj|2008t) . and other galaxies jBender et al.lll994Tl . 

h Logarithmic value of X-ray luminosity: NGC 5128 llO'Sullivan et al.ll200ll , l2003ll , and other galaxies llBeuing et al.lll999l) . 

1 Numbers of blue GCs and red GCs: NGC 5128 llWoodlev et alj|2007h . NGC 1407 l lForbes et alj|2006h . and other galaxies 
<Kim et alj|2006ll . 

j Specific frequency of GCs: M60 <Lee et all l2008bl: iForbes et al.l l2004h. NGC 5128 llHarris et all |2004| . l2006h. NGC 4636 
llKissler et al.lll994l : Ipirsch et al.ll2005l) . NGC 1407 jPerrett et all 1 199ft IForbes et alj|200rj) . and other galaxies jBrodie fc Straderl 
l2006h . 



-33- 



Table 4. Global Kinematic Properties of GCs in gEs 



Galaxy 


GC N a p 


O SIR 




SIR/ (Tp r 




(km s~*) (km s~*) 


(deg) (km s~*) 


(km s~*) 





M60 


AGC 


121 


1073 




BGC 


83 


1086 




RGC 


38 


1040 



+22 
-22 
+27 
-25 
+48 
-42 



234+^ 


2251 


-12 
■14 


141^ 




65+ ' 27 

U.DO_ 22 


223ljjj 


2181 


■16 
-23 


130±|? 




u -°°-0.30 


258±il 


2371 


18 
-19 


171+ 58 

1 1 1_46 


2401^ 


71 + - 30 



M87 


AGC 


276 


1333l| 


4141H 


681** 


172l| 


399 + l| 


n 4q+ - 12 

u -^°-o.o9 




BGC 


158 


134111 


425l| 


5911? 


1811L 7 


414l| 


0.44+8:1* 




RGC 


118 


1324+J? 


40011 


79+lf 


16511 


380+2* 


0.43+8:* 2 7 



M49 


AGC 


263 


973 




BGC 


159 


954 




RGC 


104 


999 



r20 

-18 
-32 
-27 
h31 
-25 



322; 

352 
276 



-14 

-17 
-21 
-25 
-19 
-23 



1061^ 
10211 
18211 



5411 



92 
11 



r71 

-35 
r-79 
-83 



+14 

349ll 

19 
23 



321 



278: 



O.lTl^ 7 
0.271°;?* 
0.04lg;» 



AGC 


435 


14421^ 


32311^ 


30711 


31111 


326+1*. 


o.io±S:^ 


BGC 


216 


14451^ 


3591*.? 


26llg 


6911 


3641*? 


o.i9lg:S 


RGC 


219 


143911? 


285±1| 


ol£ 


4611 


2881* 6 , 


0.1618:11 



NGC 5128 AGC 


341 


536+J 


129lf 


18411 


30111 


129lf 


0.23 


BGC 


160 


526+H 


12611 


16812? 


25li 


1291? 


0.19 


RGC 


146 


55211* 


1331^ 


19ll$° 


47lg 


1321^ 


0.36 



+0.14 
-0.13 
+0.18 
-0.29 
+0.35 
-0.44 



AGC 


238 


949+ljj 


2251*, 2 


1741S 


3711 


2261*, 2 


o.i6l°:ll 


BGC 


108 


95ll| 


2511H 


ollll 


271H 


25211^ 


o.n±8:Sg 


RGC 


130 


949+2° 


20511* 


17811 


68!£ 


203l 12 


o.33i°:? 7 



AGC 


172 


1753: 


-19 
-18 


245+H 


226+j* 


86lg 


24711 2 ; 


S c )+°'** 


BGC 


86 


1757: 


-34 
-35 


240+*! 


27111? 


8711 


23711* 


0.36+8;** 


RGC 


86 


1750: 


-29 
-22 


251+2? 


2091 2 ? 


10411 


2501*1 


o.4il8:lI 
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Table 5. Correlations between the GC systems and their Host Galaxies 



X 


Y 


a 


b 


rms 


rs 




Remarks 


logLx 


a AGC 
p,r 


97.66 


-3796 


27.77 


0.89 


0.01 


Fig 


. Q21» 


\ogL x 


BGC 
p,r 


110.28 


-4309 


32.48 


0.96 


0.00 


Fig 


. [Tj^a) 


\ogL x 


n RGC 
p,r 


84.04 


-3243 


27.13 


0.86 


0.01 


Fig 


. QUa) 


"star 


AGC 
p,r 


1.394 


-69.51 


30.71 


0.83 


0.04 


Fig 


■ Mb) 


&star 


BGC 
p,r 


1.538 


-88.77 


45.63 


0.71 


0.11 


Fig 


■ Mh) 


&star 


RGC 
p,r 


1.237 


-50.31 


23.93 


0.83 


0.04 


Fig 


■ Mh) 


My 


a AGC 

p,r 


-255.6 


-5389 


64.80 


-0.49 


0.33 


Fig 


.E2£c) 


My 


(jBGC 

p,r 


-279.0 


-5898 


71.90 


-0.43 


0.40 


Fig 


. IHc) 


My 


n RGC 
P,r 


-222.1 


-4658 


59.60 


-0.66 


0.16 


Fig 


.(He) 


Sn 


AGC 
p,r 


20.03 


150.1 


63.28 


0.57 


0.18 


Fig 


. H^a) 


e 


AGC 

p,r 


-1548 


552.2 


56.70 


-0.77 


0.07 


Fig 


• IHc) 


Nbgc/Nrgc 


uBGC l n RGC 
p,r 1 p,r 


-0.908 


2.163 


0.107 


-0.60 


0.21 


Fig 


• ESte) 


Nbgc/Nrgc 


{IR/^ C 


1.120 


-1.003 


0.091 


0.83 


0.04 


Fig 


.Mh) 


p,r 1 p,r 




-1.153 


1.556 


0.108 


-0.86 


0.01 


Fig 


. (Ha) 



a Linear fits with Y = aX+b using the bisector method (jlsobe et al.lll990l ). rs and 
cr(rs) represent the Spearman's rank coefficient and its significance, respectively. 
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Fig. 1. — Spatial distribution of NGC 4636 GCs with measured velocit i es. Ci rcles and dia- 
monds represent the blue GCs measured in Paper I and lSchuberth et al.l (120061). respectively . 
Squares and triangles indicate the red GCs measured in Paper I and lSchuberth et al.l (120061 ). 
respectively. The GCs with velocities larger and smaller than the velocity (v ga i = 928 km 
s _1 ) of the NGC 4636 nucleus are plotted by open symbols and filled symbols, respectively. 
The symbol size is proportional to the velocity deviation. The large solid-line ellipse repre - 
sents a boundary for the standard diameter D 2 $ of NGC 4636 (jde Vaucouleurs et al.lll99ll ). 
The photometric major and minor axes of NGC 4636 are represented by the dashed lines. 
The dot-dashed line circle marks a boundary (R = 260") between the inner and outer region 
used for the analysis of the radial variation of kinematic properties. 
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Fig. 2. — Radial velocities vs. position angles for all 238 GCs {top), 108 blue GCs (middle), 
and 130 red GCs (bottom). The solid curve represents the best fit rotation curve from Table 
[U and the dot-dashed horizontal line indicates the velocity of the NGC 4636 nucleus. The 
photometric minor axis of NGC 4636 is represented by the vertical arrows (G = 58° and 
238°). 
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Fig. 3. — Radial velocities vs. position angles for the GCs in the inner region (23" < R < 260" 
, left panels) and the outer region (260" < R < 926" , right panels). The best fit rotation 
curves for all the GCs, the blue GCs, and the red GCs within each region are overlaid with 
solid lines. The dot-dashed horizontal line indicates the velocity of the NGC 4636 nucleus, 
and the vertical arrows mark the position angle of photometric minor axis of NGC 4636. 
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Fig. 4. — Radial velocities vs. projected galactocentric distances for all the GCs (top), 
the blue GCs (middle), and the red GCs (bottom). Large open squares represent the mean 
radial velocities of GCs in the radial bins that are represented by long horizontal error 
bars. Their vertical error bars denote the velocity dispersions of GCs in the radial bins. 
The histograms in the right panels represent the velocity distribution of each sample. The 
dot-dashed horizontal line indicates the systemic velocity of NGC 4636. 
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Fig. 5. — Radial profiles of velocity dispersion for all the GCs (top), the blue GCs (middle), 
and the red GCs (bottom). Filled circles represent the velocity dispersion about the mean 
GC velocity (a p ) at each point, while open circles the velocity dispersion about the best fit 
rotation curve (o~ Ptr ) at the same point. The dispersion is calculated using the GCs within 
a moving radial bin (width of 2' ~ 8.52 kpc) that is represented by a horizontal error bar 
in top panel. The dotted and dot-dashed lines denote 68% and 95% confidence intervals on 
the calculation of velocity dispersion, respectively. The dashed horizontal line indicates the 
global value of velocity dispersion of GCs in each panel. 
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Fig. 6. — Projected number density profiles for all the GC candidates (top), the blue GC 
candidates (middle), and the red GC candidates (bottom) in the photometric catalog. Filled 
circles represent the GC candidates from the H ST/WFPC2 imag es, while open squares the 
GC candidates from the KPNO CT\ images (IPark et al.l l2009bl ). The solid line and the 
dashed line in each panel indicate the projected best fits using the NFW density profile and 
the Dehnen density profile, respectively, for each sample. 
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Fig. 7. — Left: -R-band surfa ce photometry of N GC 4636 de rived from KPNO im ages 



( IPark et al.ll2009bl . open circles) , iDirsch et al.l (120051 . crosses), and lPeletier et al.l (Il990l . filled 
circles). The dotted line indicates a projected best fit using eq. ([6]). Right: Three dimensional 
stellar mass density profile using the best fit model in the left panel with a constant -R-band 
mass-to-light ratio of T = 9.0 M Q L~^ Q . 
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Fig. 8. — VDPs for the stars and the GCs. Stellar VDPs are from iBender et al.l (11994 . open 
squares), and the GC VDPs are shown by filled and open circles with associated dotted and 
dot-dashed lines (from Fig. [5]). The solid line represents the stellar VDP calculated using the 
stellar mass model in Fig. [7] with a constant stellar mass-to-light ratio of To = 9.0 M Q L~j^ & 
and a stellar velocity anisotropy of (3 S = —0.3. The double-dot-dashed curve shows the 
stellar VDP with T = 8.0 M & L~^ & and (3 S = 0.0. Other lines represent the VDPs calculated 
using the same stellar mass model as above with NFW GC density profiles and velocity 
anisotropies of (3 C \ = 0.0 (short dashed line), +0.99, —99 (long dashed lines). 
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Fig. 9. — Total mass profiles of NGC 4636 derived from X-ray data. The solid line, 
short-dashed line, filled circles, open squares, and crosses represent the mass profiles de- 
rived by Loewen s tein fe Mushotzky (2003 LM0 3), B righenti fe Mat hews (1997, BM97), 



Mushotzky et al.l (119941 ). iMatsushita et al.l (119981 ). and IJohnson et al. I (120091 ). respectively. 



The dotted line represents the stellar mass profile derived using a constant i?-band mass-to 



light ratio of T = 9.0 M q Lr & from the result of IMatsushita et al.l (119981 ). The dot-dashed 



line and the long-dashed line represen t, respectively, the profil es of the total mass and the 
dark matter derived from the result of IMatsushita et al.l (119981 ). 



-45- 



500 



400 



g 300 



K 200 



100 



500 



400 



§ 300 



\A\ 200 



100 - 








I I I I I I I I I I I I I I I I I I I I I I I I I I I 



10 20 30 40 50 
R [kpc] 



10 20 30 40 50 
R [kpc] 



Fig. 10. — VDPs (a, c) and the aperture VDPs (b, d) for all the GCs. Filled circles in (a, 
c) represent the measured VDP shown in Fig. and those in (b, d) denote the measured 
aperture VDP. Associated, dotted and dot-dashed lines represent 68% and 95% confidence 
intervals on the calculation of velocity dispersion, respectively. Three smoothly curved lines, 
from a radially biased velocity anisotropy to the tangentially biased velocity anisotropy (from 
top to bottom, (3 C \= 0.99, 0, and —99), represent the VDPs calculated using the GC density 
profile of NFW (a,b) and of Dehnen (c, d). The solid lines represent the VDPs calculated 
using the total mass profile derived from lLoewenstein fc Mushotzkv I (120031). and the d ashed 
and the dot-dashed lin es represent the mass profile by Brighenti fc Mathews! ( 1997 ) and 
Matsushita et al.l (119981 ). respectively. The solid lines labeled with j3 c \ = —9.0 in (b) and (d) 
represent approximate fits for R > 20 kpc. 
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Fig. 11. — VDPs (a, c) and the aperture VDPs (b, d) for the blue and red GCs. The 



symbols are the same as in Fig. [TOJ but the dashed, the solid, anc 
represent the VDPs calculated using the mass profile derived from 


the dot-dashed lines 
Briehenti & Mathews 


(1997 


), 


Loewenstein & Mushotzkv 


(2003), and 


Matsushita et al. 


( 1998|) with the GC density 



profile of NFW, respectively. The solid lines labeled with (3 C \ = — 1.0 and -25 in (b) and (d), 
respectively, represent approximate fits for R > 20 kpc. 
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Fig. 12. — Kinematic parameters (the velocity dispersion for all the GCs, <Jp^ c , the ratio 
of the velocity dispersion between the blue GCs and red GCs, / , the ratio of the 
rotational velocity to the velocity dispersion for all the GCs, QR/a^ f ) vs. global parameter 
(X-ray luminosity, logLx, stellar velocity dispersion, a sta r, and total \^-band magnitude, 
My) for the gEs: M60 (circles), M87 (diamonds), M49 (triangles), NGC 1399 (reversed 
triangles), NGC 5128 (squares), NGC 4636 (pentagons), and NGC 1407 (stars). Open circles 
and crosses in (a), (b), and (c) represent the blue and red GCs in each gE, respectively. The 
solid lines represent the bisector linear fits. The dotted line in (b) represents one-to-one 
relation. The dashed line in (c) represents the fit without NGC 1399, which is listed in Table 
5. 
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Fig. 13. — Kinematic parameters (the velocity dispersion for all the GCs, , the ratio 

of the velocity dispersion between the blue GCs and red GCs, &p^ c l&p^ C ', the ratio of 
the rotational velocity to the velocity dispersion for all the GCs) vs. global parameters 
(specific frequency, Sn, the number ratio of the blue GCs and the red GCs, N bgc /N R gc, 
and ellipticitye) for the gEs. The symbols are the same as in Fig. [12j The solid lines 
represent the bisector linear fits. 
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Fig. 14. — The ratio of the rotational velocity to the velocity dispersion for all the GCs, 
riR/<jp^ c , vs. the ratio of the velocity dispersion between the blue GCs and red GCs, 
a p,r C /< 7 p!r C ( a )> an d the difference between the GC rotation angle and the position angle of 
the minor axis of the host galaxies, 6 — 9 minor (b) f° r the gEs. The symbols are the same 
as in Fig. [12j The solid lines represent the bisector linear fits. 
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Fig. 15. — The velocity anisotropy {(3d) vs. the velocity dispersion of all the GCs (a), blue 
GCs (b), and red GCs (c) in gEs. The velocity anisotropies only indicate representative 
positions according to the decision of each velocity ellipsoid. 
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Fig. 16. — The rotation-corrected velocity dispersions, o vx ja^ , vs. the projected galac- 
tocentric distances normalized to the effective radius, R/R e s, for all the GCs (top), blue 
GCs (middle), and red GCs (bottom) in the gEs. Open symbols indicate the dispersions in 
the inner region of each gE, while filled symbols those in the outer region. The dot-dashed 
lines represent the average of the rotation-corrected velocity dispersions for the inner region 
(R/R c s < 2) and outer region (R/R e s > 2) of each gE. 
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Fig. 17. — The ratio of the rotation amplitude to the velocity dispersion, QR/o-pf c , vs. the 
projected galactocentric distances, R/R e s, for all the GCs (top), the blue GCs (middle), and 
the red GCs (bottom) in gEs. The symbols are the same as in Fig. [161 The dot-dashed lines 
represent the average of the rotation-corrected velocity dispersions for the inner region and 
outer region of each gE. 
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Fig. 18.— The ratio of the rotation-corrected velocity dispersion between the blue 
GCs and the red GCs, &p^ C l&p^ C (a), and the ratio of the rotation amplitude to 
the rotation-corrected velocity dispersion between the blue GCs and the red GCs, 
{VtR) / a^ c / {VtR) / a^ c (b) vs. the projected galactocentric distance. The symbols are 
the same as in Fig. [161 The dot-dashed lines represent the average of the rotation-corrected 
velocity dispersions for the inner region and outer region of each gE. 



